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Effectiveness of a pre-denitrification activated sludge treatment system is governed by the kinetics of the biological 
reactions, and the hydrodynamic mixing behavior in the reactors. Achieving good mixing conditions within a reactor not 
only enhances the transfer of reactants but also ensures homogeneous environmental conditions throughout the vessel 
when required, allowing for an effective usage of the reactor’s total volume, leading to optimized, low-cost operation. In 
this work, a pre-denitrification activated sludge system performance with regards to the biological treatment of organic 
carbon and nitrogen was investigated, under two scenarios for non-ideal mixing in the anoxic reactor. The system 
performance is simulated based upon the Activated Sludge Model 1 model's biological reactions, and combining two 
non-ideal mixing two-parameter models: CSTR with bypass and dead volume, and two CSTRs with exchange. 
Performance discrepancies were then identified in the presence of non-ideal mixing. The system’s performance was 
found to be more susceptible to the presence of a dead volume/bypass scenario compared to the two CSTRs with 
material exchange scenario. Under non-ideal mixing conditions, effluent concentrations of Total Kjeldahl Nitrogen, 
organic carbon increased marginally, while effluent concentration of nitrate increased significantly. Similarly, the waste 
stream concentrations of Total Kjeldahl Nitrogen and organic carbon increased significantly as a result of an increase in 
the concentration of the heterotrophic biomass. The outcome of this study provides an insight when troubleshooting the 
operation of pre-denitrification activated sludge systems for non-ideal mixing conditions. 
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Introduction 


Activated Sludge (AS) systems are among the most widely employed biological treatment processes in wastewater treatment 
plants(Rittmann and McCarty, 2001). AS systems are engineered processes, designed and operated to sustain and facilitate the growth 
of different microorganisms. These microorganisms in their part treat wastewaters biologically by incorporating/removing the 
pollutants in their different metabolic and growth associated processes. Those pollutants are mainly comprised of dissolved and 
particulate organic carbon and nitrogen-containing compounds. The main biological processes leading to wastewater treatment taking 
place in AS systems, and the associated microorganisms responsible for them are organic carbon mineralization under aerobic and 
anoxic conditions by facultative heterotrophic bacteria, ammonia (NH3) oxidization under aerobic conditions to nitrate (NO3') by 
autotrophic nitrifying bacteria, and nitrate reduction under anoxic conditions to nitrogen gas (N2) by the aforementioned facultative 
heterotrophic bacteria (utilized instead of oxygen in the metabolic process) (Madigan and Martinko, 2006; Tchobanoglous et al., 
2007). Pre-denitrification is (PD) one approach for the operation of AS (Hellinga et al., 1999; Kim et al., 2009). The required 
environmental conditions (anoxic and aerobic) are divided between two distinctive Continuous-flow stirred-tank reactors (CSTRs), and 
the microorganisms (biomass) are circulated between them. Figure 1, illustrates the process scheme. The anoxic reactor receives both 
fresh organic carbon and ammonia containing wastewater and a nitrate-rich, biomass containing, mixed liquor recycle from the aerobic 
reactor, where a portion of organic carbon is mineralized and nitrate is reduced to nitrogen gas. Anoxic reactor effluent carries the 
remaining untreated wastewater to the aerobic reactor, where air is supplied to sustain an elevated oxygen concentration, and enhanced 
mixing conditions. The remaining organic carbon is completely mineralized and ammonia is oxidized to nitrate. The settler allows for 
biomass to be separated from the mixed liquor, and recycled back so that the system’s solids retention time is controlled independently 
from the system’s hydraulic retention time (Tchobanoglous et al., 2007). Effectiveness of a PD-AS treatment system is governed by 
the kinetics of the biological reactions, and the hydrodynamics of mixing in the reactors. The different reactions and processes taking 
place in a PD-AS system were first identified and modeled by the International Water Association (IWA) Activated Sludge Model 1 
(ASM1) (Henze et al., 2000), and since employed in a wide variety of studies (Alex et al., 2008; Hajaya and Pavlostathis, 2013; Ostace 
et al., 2011; Van Loosdrecht et al., 2015). ideal mixing conditions are commonly assumed while designing biological reactors for 
wastewater treatment processes (Vanrolleghem et al., 2003). Achieving good mixing conditions within reactors not only enhances the 
transfer of substrates (reactants) but also ensures homogeneous environmental conditions throughout the vessel, thereby, allowing for 
an effective usage of the reactor’s total volume leading to optimized and low-cost operation (Badkoubi et al., 1998; Garcia et al., 
2005). 
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Fig. 1 Typical pre-denitrification activated sludge system including anoxic and aerobic reactors and settler 


Despite the assumption of ideal mixing in biological reactors, it is likely non-ideal mixing conditions will prevail in the real situation. 
Kjellstrand (2006) studied the hydraulic behavior of a denitrifying activated sludge tank, located at the Rya Wastewater treatment plant 
(WWTP) in Goteborg Sweden, where it was found that up to 30% of the feed had notably fast hydraulic retention times, while up to 
21% of the feed had notably slow hydraulic retention times compared to the nominal hydraulic retention time. This indicated the 
presence of feed bypass and stagnation zones in the system. Sanchez and coworkers (Sanchez et al., 2016) performed the hydraulic 
characterization for the secondary treatment system in the WWTP of San Pedro del Pinatar. They reported that only 70.5% of the 
system’s volume was active, while the remaining 29.5% was actually a dead volume. Collivignarelli and coworkers (Collivignarelli et 
al., 2018) reported a dead volume fractions ranging from 15% to 25%, and bypass fractions ranging between 35% and 40%, while 
performing the experimental verification of reactor hydrodynamics for an undisclosed WWTP. Manenti and coworkers (Manentiet al., 
2018) have identified a 5% fraction of dead volume in a pilot-scale WWTP. Non-ideal mixing occurrences in AS systems were 
reported for constructed, operational wastewater treatment plants, which could be brought by poor design or equipment failure. The 
cost of diagnosing their presence in such systems could be drastically reduced if mathematical simulations are used. The purpose of 
this work is to identify the possible consequences of non-ideal mixing on the biological processes performance. Such discrepancies in 
performance can be used in identifying the presence of non-ideal mixing behavior in the AS system reactors. In this work, a pre- 
denitrification activated sludge (PD-AS) unit performance with regards to the biological treatment of organic carbon and nitrogen is 
investigated under the conditions of non-ideal mixing in the anoxic reactor. The system performance is simulated based upon the 
ASM1 model’s biological reactions (Henze et al., 2000), and combining two non-ideal mixing two-parameter models: CSTR with 
bypass and dead volume, and two CSTRs with exchange (Fogler, 1999). Performance discrepancies could then be identified by 
comparing the performance of the system under non-ideal mixing to performance of the system at ideal mixing conditions. 


1 Mathematical Models 
1.1 Biological reactions 


The ASM1 model identified the specific processes in PD-AS to be: aerobic growth of heterotrophs, anoxic growth of heterotrophs, 
aerobic growth of autotrophs, decay of heterotrophs, decay of autotrophs, ammonification of soluble organic nitrogen, hydrolysis of 
entrapped organics, and hydrolysis of entrapped organic nitrogen (Henze et al., 2000). Tables 1 detail each process mathematical 
model, where S denotes the soluble constituents’ concentrations, and X denotes particulate constituents’ concentrations, while Table 2 
lists the different variables associated with the aforementioned processes, with the kinetics of each step involved in the PD-AS system. 


1.2 Ideal completely mixed reactors model 


Figure 2 shows the PD-AS system. In this work, both the widely adapted COAST benchmark and BSM1 benchmark(Alex et al., 2008; 
Copp, 2002) are used for system size, albeit with some modifications. The system is designed for an average wastewater flow rate (Q) 
of 18400m?/day, with=300 mg COD/L of biodegradable organic carbon and ~50mg N/L of nitrogen. Volumes of the anoxic and 
aerobic reactors are 2000 m° (V7) and 4000 m° (V2), respectively, with nitrate, recycle ratio (R,) and solids recycle ratio (R2) of 3 and 1, 
respectively. The waste flow rate (Qw) is chosen for a 16-day SRT to be at 385m*/day. Air is only introduced into the aerobic reactor to 
sustain an oxygen concentration of 2.1 + 0.11 mg/l. The settler sub-model is not included in the system, and it is assumed that no 
reactions are taking place in it; instead, a particulates separation coefficient (SE) is assumed and used to calculate the concentration in 
the settler overflow (final effluent) and underflow. Particulate concentration (X) in both the settler’s overflow (Xor) and underflow 
(Xur) is determined by mass balance around the settler, by using a settling efficiency (SE), which is assumed constant (Le Moullec et 
al., 2011) and can be defined as SE=Particulate mass in the underflow/ Particulate mass in the overflow: 


15 


Jordanian Journal of Engineering and Chemical Industries (JJECI) Research Paper, Vol. 2 No.1 2019 


Xur = [SE. Q. X2. (Rz + 1)]/[R2.Q + Qw]... (1) 
Xor = [(1 — SE). Q. X2. (R2 + 1)]/[Q + Qwl... (2) 
where X2is particulate concentration (mg COD/1) in the aerobic reactor. 


Table 1 Biological processes taking place in the PD-AS, adapted from the ASM1 model (Henze et al., 2000) 


Process Model Symbol 
Ss So 
Aerobic growth of heterotrophs My ahaha BH P; 
; s OH Sno 
Anoxic growth of heterotrophs Hy 5, +K, 5, + Kou Sho + Kno NGXBH P2 
Aerobic growth of autotrophs Swn x BA P; 


Ma o ER SG + Koa 


Decay of heterotrophs byXepu P4 
Decay of autotrophs baXBa Ps 
Ammonification of soluble organic nitrogen kaSnpXBH Ps 
X;/Xpu Kon SNo Kon 
Hydrolysis of ent d ics + X P 
ydrolysis of entrapped organics "X [Xen + Ky (So + Kon Nn Sno + Kuo So + Kop! “P” 7 
Hydrolysis of entrapped organic nitrogen 7 ae Ps 
s 
Table 2 State variables in the PD-AS model and their associated rates 
Variable Definition Rate 
S? Soluble inert organics Ys, = 0 
-1 
Ss Readily biodegradable (soluble) substrate Tss = T (P, + P2) + P; 
H 
Xf Particulate inert organics Ty, = 0 
Xs Slowly biodegradable (particulate) substrate Tys = (1 — fp)(P, + Ps) — P; 
Xop Active heterotrophic biomass Tyga = Pi + Po — Py 
Xp Active autotrophic biomass Txga = Ps — Ps 
Xp" Non-biodegradable particulates Typ = fp(Py + Ps) 
1- 4.57 — 
So? Dissolved oxygen Tso == Mi P — oe P3 
Yu Va 
Sno" Nitrat Emel ree: 
Tsvo = a 
VO’ itrate sno = 2 B6y,4 j ya? 
Swat Free and ionized ammonia Tsu = (-ixe = 5) Pz — ix (P, + P2) + Pe 
A 
Swo* Soluble biodegradable organic nitrogen Tsnp = Pg — Pe 
Xyp Particulate biodegradable organic N Tyno = Cixg — İxp fp) (P4 + Ps) — Pg 
X 
Sarr? Alkalinity ua 


Tsark = Pz X 
s 


a: mg COD/L, mg of Chemical Oxygen Demand; b: mg Oz L; c: mg N/ L; d: mol/L. 


Assuming ideal CSTR behavior in the anoxic and aerobic reactors, constant liquid densities and reactor volumes, and reactions taking 
place only in the reactors, the following equations can be used to describe the dynamic behavior of soluble (S) and particulate 
constituents: 

For the anoxic reactor: 


dS,/dt = (Q/V1) ([S; + S2(R1 + R2)] — [1 + Ry + R2]S1) + fsı... (3) 
aX, /dt = (Q/V,) ([X; + X2R1 + XypR2] — [1 + Ry + R2]X1) + rı... (4) 


For the aerobic reactor: 


dS,/dt = (Q[1 + Ry + R2]/V2) (S1 — S2) + 152... (5) 
dX,/dt = (Q[1 + Ry + R2]/V2) (X1 — X2) + Tyo... (6) 


Where S; and X; are the concentrations in the anoxic reactor, S and X; are the concentrations in the feed wastewater (Table 3), and S2 
and X2 are the concentrations in the aerobic reactor, respectively. Xur is the particulate concentrations in the solids recycle (Equation 
1), and rs;«s2 and rx;ex2 are the rates in the different reactors (Table 2). Finally for oxygen, in the aerobic reactor, the following rate is 
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introduced to represent the aeration rate (FAir-in): 
Tair—in = KL (S^ — Soz) 


where Kz is the oxygen transfer coefficient (d) maintained at a level in the aerobic reactor to provide a constant oxygen concentration 
(as mentioned above)(Alex ef al., 2008), and S347 is oxygen saturation concentration =8mg/L, (at 26°C and latm) (Rittmann and 
McCarty, 2001), and So2 in the oxygen concentration in the aerobic reactor (mg/L). 


1.3 Non-ideal mixing models 


Among the PD-AS system reactors, the anoxic reactor is the most susceptible to the occurrence of non-ideal mixing. Mixing in this 
reactor is solely performed mechanically, compared to superior mixing conditions (minimizing non-ideal mixing) brought by the 
aeration process taking place in the aerobic reactor. Accordingly, the non-ideal mixing scenarios are assumed to affect the anoxic 
reactor in the PD-AS. Oxygen concentration and settler behavior are handled in the same matter as discussed above in the ideal CSTR 
based system. 


oe ae 


ep +¢- High Agitation 


Sts- --Low Agitation 


Fig. 2 Non-ideal, real mixing scenarios; A: Bypass with dead volume, and B: Rapid and slow mixing. 


1.3.1 CSTR with bypass and dead volume 


In this situation, poor mixing results in generating a dead zone within the reactor volume, where limited or no exchange of material is 
taking place between it and the remaining reactor volume. This will reduce the available volume for reactions within the reactor itself. 
In addition, poor mixing can allow for a portion of the feed to exit rapidly without being mixed inside the reactor, as if that portion is 
totally not entering the reactor, but rather bypassing it directly to the effluent (Figure 2.A). Previous scenarios are modeled as a 
combination of an ideal CSTR coupled with a dead zone. These volumes are defined as fractions of the total reactor volume (V), 
where a = volume of CSTR/V and (l-a) = volume of dead zone/V. On the other hand, flow rates headed for the reactor (Q) have a 
fraction that bypasses it while the remaining fraction enters it, where f= bypassed flow/Q and (1-f)= entering flow/Q (Fogler, 1999). 
The degree of nonideality depends on the values of a and 2. The presence of a bypass and a dead volume in real engineered biological 
reactors has been previously reported (Collivignarelli et al., 2018; Kjellstrand, 2006; Manenti et al., 2018; Sanchez et al., 2016). In 
light of the previous argument, the dynamic behavior of soluble and particulate constituents in both reactors no longer behaves as 
described by Equations 3, 4, 5, and 6. Assuming constant values for a, J (for all streams entering), liquid densities and reactor 
volumes, and reactions taking place only in the reactors, the following equations can be used to describe the dynamic behavior of 
soluble and particulate constituents: 

For the anoxic reactor: 


dS,/dt = (Q(1 — B)/aV,) ([S; + S2(R1 + R2)] — [1 + Ry + R2161) + fsı... (8) 
dX, /dt = (Q(1 — B)/aV,) ([X; + XR, + Xur R2] — [1 + Ry + R2]X1) + 11... (9) 
For the aerobic reactor: 

dS,/dt = (Q/V2)((1— B)[1 + Ry + Ro|S, + BLS; + S2(Rı + R2)] — [1 + Ry + R3]S2) + fsz... (10) 
dX,/dt = (Q/V2) ((1 — B)[1 + Ry + R2]X1 + BIX; + X2Ry + Xur R2] — [1 + Ry + Ry] X2) + 1x2... (11) 


1.3.2 Two CSTRs with material exchange 


This scenario depicts a situation where mixing apparatus are poorly situated within the reactor, resulting in rapid mixing in their 
vicinity, while the remaining region undergoes mixing, albeit at a lesser rate. This situation will influence the distribution of material 
within the reactor, affecting the different reactions rates due to an uneven distribution of substrates (Figure 2.B). In this situation, it is 
assumed that the anoxic reactor volume (V) is split into two fractions: a fraction that undergoes high agitation (a=highly agitated 
volume/V) and a fraction with low agitation ((1-a)=volume with less agitation/V). Material is exchanged between the two-volume 
fractions by a ratio of the flows entering the reactor; f= material exchange flow rate/Q. Both the inlet and outlet flow is directed to and 
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from the highly agitated portion (Fogler, 1999). As with the previous model, the degree of nonideality depends on the values of a and 
f. The two CSTR with exchange model has been previously used to describe non-ideal mixing in different biological reactors (Le 
Moullec et al., 2010; Le Moullec et al., 2011; Pereira et al., 2012). Assuming constant values for a, p (for all streams entering), liquid 
densities and reactor volumes, and reactions taking place only in the reactors, the following equations can be used to describe the 
dynamic behavior of soluble and particulate constituents: 


For the anoxic reactor highly agitated fraction (concentrations are accented with 1): 
dS*,/dt = (Q(1 + Ry + Rz)/aV;) ([S; + S2(Rı + R2)/1 + Ry + R3] — $4, + B(S? — S14)) 4151... (12) 
dX1,/dt = (Q(1 +R, + Rz)/aV,) ([X; + XR, + XypR2)/1 + Ry + Ro] — Xt, + B(X*1 — Xt) Pera (13) 


For the anoxic reactor fraction with less agitated (concentrations are accented with 2): 


dS?,/dt = (Q(1 + R, + R2) /(1 — a@)V,) (St; - SPDR eats (14) 
dX1,/dt = (Q(1 +R, + R2)B/C — a)V,) (Xt - X*,)41y2... (15) 
For the aerobic reactor: 

dS,/dt = (Q(1 + Ry + R2)/V2) (St, — S2) + Fp... (16) 
dX,/dt = (Q(1 + Ri + R3)/V2) (X11 — X2) + Tyo... (17) 


Table 4 Kinetic and stoichiometric parameters used in the simulation (Alex et al., 2008) 


Parameter Definition Value 
uu Max. specific growth rate for Heterotrophs (d°!) 4.0 
Table 3 Wastewater characteristics used in the ua Max. specific growth rate for Autotrophs (d‘') 0.5 
simulations (Vanhooren and Nguyen, 1996) Ks Half saturation constant for Heterotrophs (mg COD/L) 10 
Variable Value Kon Half saturation constant for O2 Heterotrophs (mg O,/L) 0.2 
Sif 30 Kyo Half saturation constant for Heterotrophs (mg NO3-N/L) 0.5 
Ss? 70 nc Correction for Anoxic Heterotrophic growth (-) 0.8 
Xie 52 Koa Half saturation constant for O2 Autotrophs (mg O,/L) 0.4 
Xf 200 Kya Half saturation constant for Autotrophic. (mg NH3-N/L) 1.0 
Xpyê 28 by Decay constant for Heterotrophs (d`') 0.3 
Xp’ 0.25 ba Decay constant for Autotrophs (d"!) 0.05 
Xp 5.0 ka |Ammonification rate (l.COD/mg.d) 0.05 
So? 0.25 kn Max. specific Hydrolysis rate (mg COD/mg COD biomass.d) 3.0 
Sioi 0.25 Kx Half saturation constant for Hydrolysis (mg COD/mg COD 0.1 
SNHE 30.0 biomass) 
Swot 7.0 nn Correction for Anoxic Hydrolysis (-) 0.8 
Xwnpi® 11.0 yx Heterotrophic yield coefficient (mg biomass/mg COD) 0.67 
Sarki 10.0 ya Autotrophic yield coefficient (mg biomass/mg N) 0.24 
a: mg COD/L, mg of Chemical Oxygen fp Particulate yielding biomass fraction (-) 0.08 
Demand; b: mg O,/L; c: mg N/L; d: mol/L ixg Nitrogen fraction in biomass (ng N/mg COD biomass) 0.08 
ixp Nitrogen fraction in biomass products (mg N/mg COD biomass) 0.06 


2. Computer simulations 


The performance of the PD-AS is simulated dynamically under ideal conditions (CSTR) and real conditions (CSTR with dead volume 
and two CSTRs with material exchange). The characteristics of treated wastewater are given in Table 3, the kinetic and stoichiometric 
parameters are given in Table 4. The group of ODEs representing the behavior of all constituents in the system where solved 
simultaneously in order to simulate the operation of the system for 30 days. The equations were solved using fourth-order Runge— 
Kutta procedure in MATLAB (The MathWorks Inc., Natick, MA), with a maximum time step of 1 day. The system’s performance was 
evaluated by calculating the following constituents in the effluent (Alex et al., 2008): 

A. Total Kjeldahl Nitrogen (Rittmann and McCarty, 2001) in effluent and waste (TKNz and TKNw): 


TKN; = Siu + XNp + ixe(Xpa + Xen) + ixp(Xp" + XP") (18) 
TKNw = Sin + Xnp + ixe (XBA + Xen) + ixp(Xp" + XP") (19) 
B. Total Ammonia and Nitrate concentrations (Syq-¢ and Swo-z). 
C. Biochemical oxygen demanding (Rittmann and McCarty, 2001) organics in the effluent (BOD; and BOD): 


BOD, = 0.25(S2 + XSF + (1 — f)(X8i + X84)) (20) 
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BODy = 0.25(S2 + XF + (1 — f) KET + XBF)) (21) 
D. Chemical oxygen demanding (Rittmann and McCarty, 2001) organics in effluent and waste (COD; and CODy): 

COD,g = S? + S + XP" + XEF 4+ XRF + XPE + XPF (22) 

CODy = SÈ + Sé + XPE + XU" + XPI + XPE +XRE (23) 


To study the system performance under real conditions, multiple simulations were performed while varying the values of the two 
parameters in the real reactor models in order to reflect varying degrees of nonideal flow: For the CSTR with dead volume model, a 
was changed from 1 to 0.1, while £ was changed from 0 to 0.9, and for the two CSTRs with exchange model a was changed from 0.9 


to 0.1, while J was changed from 0.01 to 0.9. The aforementioned values for a and J were chosen to depict a moderate, intermediate, 
and extreme non-ideal flow scenarios. 
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Fig. 3 PD-AS performance under ideal mixing condition; A: TKNz, B: TKNw, C: Swu-z, D: Syo-z, E: BODg and CODz, and F: BODw and CODy. 


3 Results and Discussion 
3.1 Ideal completely mixed reactors 


The PD-AS system operation was simulated for 30 days. Figure 3, shows the system performance as a function of time. Effluent 
steady-state Sy, Sno, and TKN: concentrations were 1.8, 12.5, and 3.9 mg N/L, respectively, while BOD; and COD; were 16.4 mg/L 
and 48.9 mg/L, respectively. The daily waste included 626 mg N/L of TKN , 1271 mg/L of BOD, and 8625 mg/L of COD. The 


simulation outcome was comparable to other published simulations results for similar PD-AS reactor sizes and feed wastewater 
composition (Alex et al., 2008). 


3.2 CSTR with bypass and dead volume 
Figure 4 shows the effect of varying the non-ideal mixing parameters on the concentration of nitrogen compounds in the system. As 
can be seen, the system’s effluent TKN concentration varied marginally with increasing mixing non-idealities (Figure 4A) (ie. 


increased bypass fraction and dead volume fraction). This can be attributed to the fact that effluent TKN is mainly soluble, and 
removed in the aerobic reactor, where mixing performance was assumed ideal. 
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Fig. 4 Simulated PD-AS nitrogen related performance under CSTR with bypass/dead volume non-ideal mixing effects at 0.1 < æ < 1.0 and0 < f < 0.9; A: 
TKNr, B: TKNy, C: SNH-E, and D: SNo-E- 


The same trend was noticed for the total ammonia in the effluent, which remained relatively constant (Figure 4C). TKNw increased 
with increasing mixing non-idealities, as can be seen in Figure 4B. This is attributed to the increase in the heterotrophic biomass in the 
system (from 2570 to 2980 mg COD/L) as a result of the increase in substrate availability in the aerobic reactor brought by the increase 
in bypass ratio. The TKNy increase may affect anaerobic treatment for the wasted biomass (Chen et al., 2008; Tezel et al., 2014), due 
to the elevated levels of nitrogen. On the other hand, the effluent nitrate concentration increased drastically as the system was shifted 
towards non-ideal mixing conditions. As seen in Fig. 4C, Syo.¢ exceeded 20 mg N/L at a < 0.4 and p < 0.4. The decrease in a 
(increase in dead volume fraction) had a larger effect on Sno- compared to J. Similar performance short comes were observed by 
Manenti et al. (Manenti et al., 2018), when the actual retention time failed in assuring an acceptable treatment level below a minimum 
threshold value. Collivignarelli (Collivignarelli et al., 2018) reported that reducing the dead volume fraction will result in an 
enhancement of nitrate removal, while Kjellstr and co-researchers (Kjellstrand et al., 2005) reported that the presence of a dead volume 
in the reactor will reduce the denitrifying capacity due to reduced active volume, and high nitrate concentration in the effluent can 
appear due to the short-circuiting stream (bypass). Figure 5 shows the effect of varying the non-ideal mixing parameters on the 
concentration of BOD and COD in the system. As seen in Figure 5A and Figure 5B, increasing mixing non-idealities had a slight 
effect on BOD; and COD;. However, as with Syo-z, decreasing a had the main effect. This can be attributed to the reduction in the 
anoxic reactor's active volume size, leading to a lesser conversion due to smaller real retention times (Fogler, 1999). This connects to 
the increase in Syo.z as seen in Figure 4D since (as mentioned in section 1) nitrite is utilized exclusively by the heterotrophic biomass 
organic carbon mineralization in the anoxic reactor. 
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Fig. 5 Simulated PD-ASorganic carbon-related performance under CSTR with bypass/dead volume non-ideal mixing effects at 0.1 < a < 1.0 and0 < f < 0.9; 
A: BODr, B: CODr, C: BODy, and D: CODy. 


As seen in Fig. 5C and Fig. 5D, both of BODw and CODw concentrations increased in the waste stream as a direct response to 
increased biomass in the system, which was brought by the greater than before availability of substrate in the aerobic reactor due to 
increased bypass fraction and decreased active volume fraction in the anoxic reactor. Collivignarelli (Collivignarelli et al., 2018) 
reported that reducing the dead volume fraction will enhance COD and BOD. The increase in the concentrations of BODw and CODw 
in the waste stream, coupled with the increase of TKNw, will amplify the load on waste disposal facilities in the wastewater treatment 
plant. Similar to what has been observed with Syo.~ above, the active volume fraction (a) had a larger effect on the system’s BOD and 
COD concentrations when compared to the bypass fraction. It must be added that, apart from nitrate concentration, the resulting 
changes in the PD-AS’s performance were related to particulate constituents in the system. However, their concentration in the system 
is directly correlated with the settler’s performance (i.e. SE value). The SE value was assumed constant to eliminate the effect of the 
settler’s performance on the predicted outcome. 


3.3 Two CSTRs with material exchange 


As mentioned in section 2.3.2, this model describes a scenario where the mixing apparatus in the anoxic reactor are poorly positioned, 
resulting in rapid mixing in their vicinity, while the remaining region undergoes mixing, albeit at a lesser rate. Carried out simulations 
for the system’s operation at the current non-ideal mixing scenario showed that, to some extent, its performance was unaffected. This 
was in a specific range of the mixing model a and p values: 0.4 < a < 0.9 and0.1 < p < 0.9. Table 5 shows the unaffected system’s 
performance. Figure 6 and Figure 7 shows the performance of the system beyond the abovementioned range for the mixing model 
parameters values (i.e. second range at 0.1 < a < 0.4 and0.01 < f} < 0.1). It is clear that the system’s effluent TKN concentration 
slightly (Figure 6A) varied with increasing mixing non-idealities within the second range (i.e. reducing the rapid agitation zone and 
material exchange fractions). This took place because TKN is mainly removed in the aerobic reactor, which its mixing performance 
was assumed ideal. The relatively constant total ammonia in the effluent (Figure 6C) suggests that that ammonia removal is 
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unaffected. 


Table 5 Effluent constituents of the PD-AS system with and without ideal mixing conditions 


Effluent Constituents Ideal mixing CSTRs with material exchange 
TKN: 3.9 mg N/L 3.9+0.09°mg N/L 

TKNw 626.5 mg N/L 628.5+4.8mg N/L 

SNH-E 1.9 mg N/L 1.8+0.05mg N/L 

SNo-E 12.6 mg N/L 14.1+1.2mg N/L 

BODe 16.2 mg BOD/L 16.7+0.9 mg BOD/L 

CODe 48.9 mg COD/L 49.1+0.2 mg COD/L 

BODw 1271.9 mg BOD/L 1278.8+15.2 mg BOD/L 
CODw 8625.6 mg COD/L 8652.3461.5 mg COD/L 


a: Average + Standard Deviation. 
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Fig. 6 Simulated PD-AS nitrogen related performance under the 2 CSTRs with material exchange non-ideal mixing effects at 0.1 < a < 0.4 and0.01 < B < 
0.1; A: TKNez, B: TKNw, C: SNH-E, and D: SNO-E- 


However, the TKNw (which is mainly comprised from particulates) increased with increasing mixing non-idealities (Figure 6B), but 
only at extreme non-ideal mixing condition, when the anoxic reactor volume was mostly affected by low mixing due to very slow 
agitation (a < 0.3 and $ < 0.01). Operation within this range also increased the heterotrophic biomass in the system (from 2570 to 
2708+136 mg COD/L), due to ineffective utilization of readily biodegradable organic carbon in the anoxic reactor (Ss+Xs increased in 
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the anoxic reactor effluent from 82.2 to 117412mg COD/L), and were utilized by the aerobic reactor. Similar to the dead 
volume/bypass model, the increase TKNw concentration and biomass concentration in the waste stream may affect anaerobic treatment 
processes for the wasted biomass. The effluent nitrate concentration increased drastically as the system was shifted towards extreme 
non-ideal mixing conditions. As seen in Fig. 6C, Syo.e exceeded 20mg N/L at æ < 0.2 and £ < 0.01. This behavior mimicked the 
result seen in the dead volume/bypass model, however at more extreme non-ideal mixing conditions in this model. The increase is 
probably a result of the anoxic reactor's active volume becoming extremely low (i.e. fast residence times), not permitting any effective 
reactions. This can be connected to the increase of readily biodegradable organic carbon in the anoxic reactor (see before), being not 
utilized fully under these conditions. Results shown in Fig. 6 showed that the system’s performance was more affected by the high 
agitation fraction (a). Fig. 7 shows the effect of varying the nonideal mixing parameters on the concentration of BOD and COD in the 
system. 
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Fig. 7 Simulated PD-AS nitrogen related performance under the 2 CSTRs with material exchange non-ideal mixing effects at 0.1 < a < 0.4 and0.01 < B < 
0.1; A: TKNeg, B: TKNy, C: SNH-E, and D: SNO-E- 


As seen in Fig. 7.A and Fig. 7.B, only extreme levels of mixing non-idealities had an effect on BODz and CODz. The increase is 
probably related to the increased heterotrophic biomass in the system (from 2570 to 27084136 mg COD/L), due to ineffective 
utilization of readily biodegradable organic carbon in the anoxic reactor. As seen in Fig. 7.C and Fig. 7.D, both of BODw and CODy 
concentrations increased in the waste stream as a direct response to increased biomass in the system, which was brought by the greater 
than before availability of substrate in the aerobic reactor due to reduced active anoxic reactor volume. The increase in the 
concentrations of BODw and CODy in the waste stream will amplify the load on waste disposal facilities in the wastewater treatment 
plant. The aforementioned results showed that the system’s performance was more affected by the high agitation fraction (a). Results 
shown in Figs. 6 and7 showed that the non-ideal mixing scenario depicted by the two CSTRs with the material exchange had minor 
effects on the performance of the PD-AS system. Even so, these effects occurred at extreme non-ideal situations. However, this 
stipulation could be related to the specific operating conditions of the PD-AS under investigation. Moreover, previously published 
work suggests that increasing the number of CSTRs (i.e. compartments) in this model could provide a more representative depiction of 
the mixing non-idealities in AS systems (Liotta et al., 2014). 
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Conclusions 


The system’s performance was largely affected by the presence of a dead volume/bypass scenario compared to the two mixing zones 
with a material exchange scenario. Under non-ideal mixing conditions, effluent concentrations of Total Kjeldahl Nitrogen, organic 
carbon increased marginally. On the other hand, the effluent concentration of nitrate increased significantly under the same conditions, 
indicating a reduction in nitrate removal efficiency. The heterotrophic biomass concentration in the system increased under non-ideal 
mixing conditions. This by its part resulted in an increase in the waste stream concentrations of Total Kjeldahl Nitrogen and organic 
carbon. This study provides an insight into the behavior of pre-denitrification activated sludge systems when affected by non-ideal 
mixing conditions. Additionally, it identified some performance discrepancies that could-if found in a real system-indicate the presence 
of a dead volume/bypass or two mixing zones mixing non-idealities in the anoxic reactor. 


Nomenclature 


Acronyms 

AS =Activated Sludge 

PD =Pre-Denitrification 

CSTR =Continuous stirred-tank reactor 

SE =particulates separation coefficient 

COD =Chemical Oxygen Demand 

BOD =Biochemical Oxygen Demand 

TKN =Total Kjeldahl Nitrogen 

E =Effluent 

W =Waste 

UF =Underflow 

OF =Overflow 

Symbols 

S =Dissolved material concentration 

Si =Dissolved material concentration in the feed 
Sı =Soluble inert organics concentration 

Ss =Readily biodegradable (soluble) substrate concentration 
X =particulate material concentration 

Xi =Particulate material concentration in the feed 
Xı =Particulate inert organics concentration 

Xs =Slowly biodegradable (particulate) concentration 
Xpy =Active heterotrophic biomass concentration 
XBA =Active autotrophic biomass concentration 
Xp =Non-biodegradable particulates from cell decay concentration 
So =Dissolved oxygen concentration 

Sot =Saturation dissolved oxygen concentration 
SNo =Nitrate concentration 

SNH =Free and ionized ammonia concentration 

SND =Soluble biodegradable organic nitrogen concentration 
Xnp =Particulate biodegradable organic nitrogen concentration 
SALK =Alkalinity 

Kı =Oxygen transfer coefficient 

Q =Feed wastewater flow rate 

Ow =Waste flow rate 

R; =Nitrate recycle ratio 

Rə =Solids recycle ratio 

Vı =Anoxic reactor volume 

V2 =Aerobic reactor volume 

r =Reaction rate 

Ks =Half saturation constant for Heterotrophs 
Kon =Half saturation constant for O2 Heterotrophs 
Kyo =Half saturation constant for Heterotrophs 
Koa =Half saturation constant for O2 Autotrophs 
Kyu =Half saturation constant for Autotrophic 

bu =Decay constant for Heterotrophs 

ba =Decay constant for Autotrophs 

Ky =Ammonification rate 

kn =Max. specific Hydrolysis rate 

Kx =Half saturation constant for Hydrolysis 

YH =Heterotrophic yield coefficient 

YA =Autotrophic yield coefficient 

Tp =Particulate yielding biomass fraction 

ixB =Nitrogen fraction in biomass 

ixp =Nitrogen fraction in biomass products 


[mg COD/L] 
[mg COD/L] 
[mg COD/L] 
[mg COD/1] 
[mg COD/L] 
[mg COD/L] 
[mg COD/L] 
[mg COD/L] 
[mg COD/L] 
[mg COD/L] 
[mg COD/L] 
[mg OXL] 
[mg OXL] 
[mg N/L] 
[mg N/L] 
[mg N/L] 
[mg N/L] 
[mol/L] 


[mg/L.d] 

[mg COD/L] 

[mg O2/L] 

[mg NO3-N/L] 

[mg O2/L] 

[mg NH3-N/L] 

[d] 

[d"] 

[L.mg COD/mg.d] 

[mg COD/mg COD.d] 
[mg COD/mg COD] 

[mg biomass/mg COD] 
[mg biomass/mg N] 

[-] 

[mg N/mg COD biomass] 
[mg N/mg COD biomass] 
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Greek letters 


a =Non-ideal mixing parameter 1 [-] 
B =Non-ideal mixing parameter 2 [-] 
Hu =Max. specific growth rate for Heterotrophs [d] 
HA =Max. specific growth rate for Autotrophs [d] 
NG =Correction for Anoxic Heterotrophic growth [-] 
Nh =Correction for Anoxic Hydrolysis [-] 
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